Single electron spin and its coherence in Si quantum computer architecture 
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The possibility of performing single spin measurements in Si-based quantum computers through 
electric field control of electrons bound to double donors near a barrier interface is assessed. We find 
that both the required electric fields and the tunneling times involved are probably too large for 
practical implementations. On the other hand, operations with double donors in their first excited 
state require smaller fields and faster tunneling times, and are therefore suitable for spin-to-charge 
conversion measurements. We also propose a measurement scheme that would render statistical 
(ensemble) estimates of the spin coherence at the Si/Si02 interface. 
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Among the several operations required for a spin-based 
quantum computer (QC), single spin rotations and mea- 
surement are probably the hardest ones to achieve within 
the currently available technology. Not only is the elec- 
tron spin extremely weak to be detected operation 
and measurement times must be fast enough compared 
to electronic spin dephasing times. Single electron spin 
control has very recently been demonstrated in a double 
quantum dot (QD) configuration in a GaAs heterostruc- 
ture by Koppens et al. [2( • In the occupation number rep- 
resentation for the individual dots, the involved states are 
(1,1), one electron in each dot, and (0,2), one of the dots 
empty and the other doubly-occupied. Due to Pauli prin- 
ciple, state (0,2) is only accessible if the electrons form 
a spin singlet, which allows inferring spin states from 
charge transport measurements |2|. Although a similar 
experiment has been suggested using double donors in 
Si 0] and other approaches have been attempted 0] for 
single donor-based QC in Si 0, progress in Si systems has 
been much slower than in the corresponding GaAs quan- 
tum dot systems, primarily because precise quantitative 
practical schemes for spin-to-charge conversion have not 
been theoretically proposed in Si QC architectures. The 
proposal by Kane et al. [| has many of the basic ingredi- 
ents that led to the recent successful results in GaAs QDs: 
The double dot potential would correspond to a double 
well structure formed by a double donor in Si near the in- 
terface with a barrier (e.g. Si/Si02) under a uniform elec- 
tric field applied perpendicular to the interface. The (1,1) 
state would correspond to one electron at the characteris- 
tic triangular shaped interface-plus-electric-field well and 
one at the double donor Coulomb well, while the (0,2) 
state would correspond to both electrons bound to the 
(neutral) double donor. 

We explore here the use of double donors in Si (solid- 
state analogues of the He atom) to perform one-qubit 
operations similar to those reported in Koppens's exper- 
iments. For the measurement of a single spin state, we 



quantitatively estimate the required electric fields for the 
first ionization of Te in Si, assessing the practical imple- 
mentation of the scheme in Ref. y. We also present a 
related proposal for the measurement of spin coherence 
times in Si near a Si02 interface. This is a key parameter 
if qubit measurement in Si-based QC is to be performed 
at the interface ||. The idea proposed here is similar 
to the one that has allowed the measurement of spin de- 
phasing times in double quantum dots with two electrons 
in GaAs In this case we show that S double donors 
would be the most appropriate due to the negligible spin- 
orbit coupling. 

Double donors (S, Se, T) in Si are substitutional deep 
centers whose electrons' binding energies (summarized in 
Table [J are typically one order of magnitude larger than 
for single donors (P, As, Sb). The substitutional dou- 
ble donor ground state is a spin singlet analogous to the 
He atom ground state. For comparison, we recall that 
the two-electrons binding energy in neutral He is 5.7 Ry, 
while the one-electron binding energy in He + is 4 Ry, thus 
the first ionization energy in He is 1.7 Ry 0. The partic- 
ular band structure of Si and the local tetrahedral sym- 
metry of the potential leads to an orbital ground state for 
substitutional donors of Ai symmetry. The first excited 
state for double donors is such that one of the electrons is 
in the ground state and the other on the next-in-energy 
ls(T"2) orbital state [g. This outer electron experiences 
a closely hydrogenic effective potential because the dou- 
ble donor core charge is virtually screened by the inner 
ls(Ai) electron which has a very small effective Bohr 
radius. The corresponding spin states are always singlet 
for the ground state, while the first excited configuration 
can be in a triplet or a singlet, with the triplet lower in 
energy due to Hund's rule. The binding energies of the 
ground and the first excited states of S, Se, and Te are 
illustrated in Fig. ^ 

We consider a double donor (S, Se or Te) on its ground 
state in Si, a distance d away from an interface with Si02- 
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Donor 


E D (meV) 


E D+ (meV) 


I s ionization (meV) 


Te 


609.6 


410.8 


198.8 


Se 


899.9 


593.3 


306.6 


S 


931.8 


613.5 


318.3 



TABLE I: Energy of neutral (Ed) and singly ionized (E D +) 
double donors. The first ionization energy is E d +-Ed. Data 
taken from Ref. lid 
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FIG. 1: (Color online) First ionization energies of the Is 
levels of double donors (Te, Se, S) in Si. The ground state 
is the singlet ls( 1 Ai)= [ls(Ai), ls(Ai)] and the first excited 
states are ls( ! T2) = [ls(Ai), ls(T2)]. (i ~ 1) corresponds to a 
singlet while (i — 3) to a triplet. Data taken from Refs. HfllTlI 



When a characteristic electric field F c perpendicular to 
the interface is applied, one of the electrons may tun- 
nel towards the interface well (where it still feels the 
Coulomb attraction of the donor which keeps it from 
spreading in a 2-dimensional electron gas [9J). Small 
changes of the electric field around F c cause the electron 
to move between the donor and the interface. This 'shut- 
tling' is only allowed as long as the two electrons are in a 
spin-singlet state while if they are in a spin-triplet state 
the electron cannot tunnel back to the donor due to Pauli 
principle. The electron motion can be detected by single- 
electron transistors (SETs) located on the device surface, 
above the oxide layer [3J . As a result, the measurement of 
charge directly leads to information about the spin state 
of the electron (spin-to-charge conversion). 

As shown in Table HJ Te is the shallowest of the double 
donors in Si, with a total (two electrons) binding energy 
i?To = 609.6 meV. For this reason it is the most apro- 
priate for electric field driven ionization experiments in 
Si, as it requires the smallest electric field F c To to ionize. 
This field has been estimated, within a single electron ap- 
proach, in Ref. 3). Here we consider the full two-particle 
Hamiltonian, which includes the electron-electron inter- 
action and the image charges, to calculate this ionization 
field as a function of d. 

We solve the 2-electron double donor problem within 
several simplifying approximations, namely single- valley 
effective mass approach, assuming isotropic effective 
mass m*. The validity of these approximations for ob- 
taining ionization fields is discussed below. The full two 
particle Hamiltonian in atomic units (a* = h 2 esi/m*e 2 



and Ry* = m*e /2h e Si ) for a double donor in Si close to 
a Si/SiC>2 (0,0,1) interface in the presence of an electric 
field is 
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with kinetic energy Tj = —d 2 /dxf - d 2 /dyf - d 2 /dzf 
where i = 1, 2, k = 3.89 x 10~ 7 e| l (m/m*) 2 cm/kV, and 
the electric field F is given in kV/cm. We estimate the 
isotropic radius ax c and mass m* from the experimental 
value of the binding energy of Te in the ground state, 
Exe — —609.6 meV, using the calculated expression for 
He, E Te = 2/a| c - 8/a Te + 5/4a Te = -5.6953i?y* which 
is minimized by a Tc = 0.59a* Q. This gives Ry* = 107 
meV, m* ~ m, and a* — 0.6 nm. The image terms are 
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the interaction of the electrons with the double donor 
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the interaction of the electrons with their own 
images and the other electron's image. Q — 
(esi0 2 - £Si)/(esi0 2 + S3i), where e Si = 11.4 and e Si o 2 = 
3.8. 

We solve the Hamiltonian in the non-orthogonal basis 
defined by the initial state = (0, 2) with the two 
electrons on the donor ground state [ls(Ai),ls(Ai)]: 



^ = ^(1)^(2) , 



(4) 



and the final singly ionized donor state and one electron 
at the interface ^>b = (1,1) 



1 



v/2(l + SJd) 



hMl)M2) +^(2)^(1)] , (5) 



where i = 1, 2, i/jd (i) ^ e ri ^ aTa , tpi(i) oc (zi + 
^2 e - Q ( 2l +d)/2 xe -/3 2 P f/2 j wdSlD = tyjjtyj) is the over- 
lap between the single electron states. The variational 
parameters a^ 0l a and minimize the expectation value 
of the energy and lead to approximate expressions for the 
one-electron ground state ipo and ipj at the donor well 
and at the interface well respectively. 

Note that the Bohr radius of the double donor ground 
state (axo = 0.354 nm) is one order of magnitude smaller 
than the one for the single donor (dp = 2.365 nm) due 
to the larger binding energy. This fact allows two use- 
ful simplifications to be made. First, the overlap Sid is 
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FIG. 2: (Color online) (a) Critical electric field F? e (cir- 
cles) required to take one of the two donor electrons to the 
interface. _F c Te depends strongly on the distance d of the 
donor from the interface. We compare it to the solution to 
the single donor problem (P) in Ref. (solid line) where an 
anisotropic donor wave-function was used. The result for an 
isotropic donor wave-function for P (dot-dashed line) with 
m* = 0.29819m gives a very similar result to the more real- 
istic anisotropic one. (b) Tunneling time required to ionize 
an electron on the first excited state of a double donor. This 
time is the same as required to ionize the single donor P [Sj. 



very small, therefore the exchange part of the electron- 
electron interaction can be neglected. Second, in the 
limit in which d >> a^ c the electron at the interface 
sees a nuclear image with charge +1 (H-like atom) rather 
than a +2 donor plus an electron (first ionized He-like 
atom). Therefore, all the image terms (Eqs. [21 and |3J) in 
the Hamiltonian are reduced to 
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and the variational parameters for ipi(i) are the same as 
for the single donor problem ||. 

We calculate F^f e , the electric field required to take 
one donor electron to the interface, from the condition 
($A\H\y A ) = {^b\H\^ b )- In Fig. Ufa) we show F c Te 
versus d. The results are similar to the single electron 
approach Q : when the double donor is located a dis- 
tance d = 20 nm from the Si0 2 interface, F c Tc = 190 
kV/cm. These values are very large and would prob- 
ably cause electrical breakdown in the nanostructures. 
The lower curves in Fig. Efa) correspond to the single 
donor P results || and are shown for comparison. The 
curves correspond to isotropic and anisotropic d wave- 
functions, respectively, and are very similar, validating 
the isotropic approximation adopted here for the double 
donor problem. In Fig. b) we reproduce the tunneling 
times Tp required to ionize P 0- T is proportional to 
the overlap between states A and B, which is orders of 
magnitude smaller in the double donor problem than in 
the single donor problem. Therefore, we expect tt c to 



be orders of magnitude larger than rp and hence prob- 
ably of the same order of the spin decoherence time in 
Si (~ 1 ms in the bulk). Consequently, double donors 
in the ground state are of no practical use to single spin 
measurements. 

This is not the case if the double donor is in its first ex- 
cited [ls(Ai),ls(T 2 )] state, which would be more directly 
accessible for electric field driven (1,1) <-> (0,2) manip- 
ulations. In contrast to the large binding energy and 
small Bohr radius of the deep centers ground state, the 
outer electron in the first excited states of neutral double 
donors has similar properties (binding energy ~ 32 meV 
and Bohr radius ~ 2 nm) || than an electron on a single 
donor ground state Therefore, F c required to singly 
ionize [ls(Ai),ls(T 2 )] and the tunneling times involved 
in this process are similar to the ones calculated for P, 
shown in Fig. [5] and, therefore, experimentaly meaning- 
ful |9j . We propose a scheme to measure spin dephasing 
times T 2 * which involves double donors in their first ex- 
cited state. As opposed to previous works which have 
measured spin dephasing times of electrons at donors in 
bulk Si 01 an d at donors 50 nm from the interface [l^ . 
our proposed experiment would measure the spin dephas- 
ing time at the interface. 

In Fig. ^ we show the first ionization energies for the 
ground state (singlet) and the first excited state (singlet 
or triplet) of the neutral double donors. The first excited 
states have been observed by looking at the ab sorp tion 
spectra from the ground state [ls(Ai),ls(Ai)] [ll|, [ToJ. 
Transitions to pure spin-triplet states are forbidden by 
selection rules however triplets have been observed in Te 
and Se 0, : Fhe spin is not a good quantum number 
due to spin-orbit interaction. On the other hand, only 
one line, due to the spin-singlet state, has been observed 
in S, due to the much smaller spin-orbit interaction for 
the lighter donor. The spin-triplet-state lifetimes in To 
and Se are of the order of tens of picoseconds as esti- 
mated from the linewidths (which could be limited by 
experimental precision) [l5T |. not much longer than the 
spin-singlet-state lifetimes. The spin-singlet lifetime in 
S is expected to be of the same order as for Te and Se, 
while the spin-triplet lifetime must be much longer due to 
the relatively small spin-orbit interaction . Moreover, 
it has been recently reported yjj that the linewidths of 
absorption transitions in P-doped Si can be decreased 
(and, therefore, the lifetimes increased) by using isotopi- 
cally purified Si, while it is expected that the linewidths 
of deeper impurities (like S) could be even more dra- 
matically reduced due to the smaller Bohr radius [l8j |. 
Therefore, in S the first excited triplet state lifetime is 
expected to be much longer than the first excited singlet 
state lifetime. In what follows, we consider Si doped with 
S. 

For S-doped Si, it has been observed experimentally 
that about 50% of the donors remain singly ionized [l0| 
(their only electron is in the inner ls(Ai) state) while 
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FIG. 3: (Color online) Double donors in their first excited 
state can be used to measure the spin-dephasing times in Si 
close to an interface. See text for details. 



times for the electron shuttling [as given in Fig.^b)], and 
the frequency of the electric field dithering (with corre- 
sponding times > 1 ns) . The proposed experiment would 
be successful in measuring spin decoherence at the inter- 
face if the spin-triplet lifetime is longer than the other 
two timescales. In this case, such measurements would 
also give estimates for the donor first excited state spin 
triplet lifetime. It is reasonable to assume that the spin 
coherence time at the interface can be made much longer 
than any of the other timescales by improving the in- 
terface quality and by isotopically purifying the Si host. 
Note that the already long bulk value ~ 1 ms can be 
dramatically increased rather easily to 100 ms or longer 
by isotopic purification 
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the others remain neutral and inert for the proposal pre- 
sented here. We could use such singly ionized donors 
to measure the spin dephasing time by manipulating 
its outer electron between the interface and the sulphur 
ls(T2) orbital (under a magnetic field < 1 T to define a 
quantization axis) in the way illustrated in Fig. 02 (1) an 
electric field slightly larger than the F c required to ionize 
the donor in its first excited state is applied; (2) elec- 
trons are injected at the interface. The ionized donors 
will bind the new electrons at the interface || leading to 
a (1,1) configuration; (3) the electric field is decreased so 
that, at F c , the electron at the interface goes to the outer 
ls(T2) state at the donor [configuration (0,2)] with typi- 
cal tunneling times as shown in Fig.[2b). (4a) If the two- 
electron spin state is a singlet, or contains a significant 
singlet component, the orbital state will rapidly decay to 
the ground state [ls(A 1 ),ls(A 1 )] ) (4b) while if they form 
a triplet, selection rules imply a much longer lived state. 
The electric field is increased again so that the outer elec- 
tron at the donor can go back to the interface. Only those 
electrons whose spins form a triplet with the inner elec- 
tron at the donor will remain in the excited state long 
enough to be able to shuttle [(1,1) <-> (0,2)] when the elec- 
tric field dithers around F c . The inner electron remains 
bound to the donor because the field required to doubly 
ionize the donor is much larger (>> F^ e ). Detection of 
charge at the interface will give a decreasing population 
of electrons from which we can extract the spin dephasing 
time of electrons at the Si/Si02 interface. 

There are various timescales involved in this problem: 
the unknown spin-triplet lifetime of the first excited state 
in S (which we have argued should be orders of magnitude 
longer than the spin-singlet ~ ps lifetime) , the tunneling 
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